Electron spin manipulation and resonator readout in a double quantum dot 

nano-electromechanical system 
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Magnetically coupling a nano-mechanical resonator to a double quantum dot confining two electrons 
can enable the manipulation of a single electron spin and the readout of the resonator's natural 
frequency. When the Larmor frequency matches the resonator frequency, the electron spin in one of 
the dots can be selectively flipped by the magnetised resonator. By simultaneously measuring the 
charge state of the two-electron double quantum dots, this transition can be detected thus enabling 
the natural frequency of the mechanical resonator to be determined. 



I. INTRODUCTION 

Nano-electromechanical systems (NEMSs) enable the 
coupling of electronic and mechanical degrees of 
freedomii^ii^ and provide an arena in which quantum 
mechanical behaviour can be identified on macroscopic 
scaleaii^i^i^i^i^. A key requirement for the realization of 
quantum electromechanical systems is the efficient trans- 
duction of the miniscule displacements of the mechanical 
element. Meanwhile theoretical studies have proposed 
coupling NEMSs to quantum two-level systems realised 
in a variety of physical systems^i^i^iiS. The coupled sys- 
tem can potentially enable the implementation of proto- 
cols similar to those realised in cavity QED and on-chip 
circuit QED with superconducting qubits and transmis- 
sion linesii. Furthermore, manipulating a two-level sys- 
tem via NEMSs has the advantage that a change in its 
state can be used to distinguish between resonance and 
off resonance conditions of the mechanical element with 
great sensitivity. 

In this work we present a proposal for coupling a NEM 
system to the spin of an electron in a double quantum dot 
(QD). Quantum dots are artificial structures which en- 
able the confinement of a few or even one electro n^^i^^'^'* . 
In this regime, the spin of the individual electron forms 
a natural two-level system and provides great potential 



for the realization of a solid-state qubit^^. However, for 
a QD based spin qubit to be realized requires precise 
control and manipulation of the electron spins. Recently, 
coherent spin manipulation in a double quantum dot was 
demonstrated by rapid electrical control of the exchange 
interaction^^ and by generating a local electron spin res- 
onance (ESR) magnetic fieldii. 

Here we describe how to exploit the versatility of a dou- 
ble quantum dot in order to realize an electron spin in a 
well defined state which is then coupled and controUably 
manipulated by the nano-mechanical resonator. Previ- 
ously, coupling spin to a mechanical resonator has been 
both theoretically and experimentally considere d^?' 
but the motivation of these studies was spin detection 
rather than spin manipulation. In the present scheme, 
the magnetised resonator selectively flips one of the two 
confined electron spins in the double dots (DDs) where 
the ESR field is produced by the oscillations of the mag- 
netised resonator. The corresponding change in the two 
electron spin state is then detected by spin blockade 
charge sensing at the DDs. Concurrently, this also en- 
ables the resonator's natural frequency to be measured. 
As demonstrated in refsJ^iil', spin-blockade readout does 
not rely on large Zeeman fields as in the case of a sin- 
gle quantum doi^ which is an important criterion for 
frequency matching between the electron spin and the 
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FIG. 1; (colour on-line) A schematic of the proposed cou- 
pled mechanical resonator double QD system. The mechan- 
ical oscillator is a suspended doubly clamped beam located 
a distance, do, from the mesa on which the double QDs are 
defined via electrostatic gates (orange). The beam is actu- 
ated to resonance via the piezoelectric effect between the top 
gate, located at where the beam is clamped (orange), and 
the shallow buried 2DEG. A micro-magnet (yellow) of mag- 
netisation, M, is incorporated onto the beam oscillator and 
moves in unison with the out of plane beam motion. The me- 
chanical oscillator couples to the confined electron spins (red 
and green arrows) in the doubles QDs via the micro-magnet's 
stray magnetic field where the degree of coupling is controlled 
by the tuning field, Bo- The quantum point contact current, 
Iqpc, is used to measure the charge state of the QDs. 

mechanical resonator. Furthermore, the effects of the 
hyperfine interaction on the electron spin states and the 
corresponding line width of the proposed measurement is 
also theoretically investigated. 



II. DEVICE GEOMETRY AND PARAMETERS 




0.1 0.2 0.3 

d(nm) 



FIG. 2: (a) The profile of the in-plane (Zeeman) magnetic 
field, Bz, near the nano-mechanical resonator with external 
tuning field. Bo = T where di and d2 mark the positions 
of the left and right quantum dots respectively. The inset 
shows the cross-section of the proposed system. The position 
of the micro-magnet is in the same plane as the 2DEG in the 
mesa on which the double QDs are located. This alignment 
is achieved by piezoelectrically flexing the beam with the ap- 
plication of DC voltage to the gate defined on the mechanical 
oscillator and is discussed in detail below, (b) The profile 
of the transverse magnetic field, B™'"', whilst the resonator 
is at maximum deflection (5 nm). The inset shows the time 
dependence of the ESR field at the positions of the quantum 
dots whilst the mechanical element is at resonance. 



A schematic of the proposed coupled nano-mechanical 
resonator DD system is shovifn in fig. [T] A GaAs / AlGaAs 
nano-mechanical beam of length 2 /im, width 500 nm 
and thickness 200 nm is located a distance, do, from the 
mesa on which the DDs reside. The natural frequency for 
a resonator of these dimensions is expected to be more 
than 200 MHz and peak displacement, at the centre of 
the beam, before the onset of non-linearity of ~ 5 nm 
assuming that the resonator can be driven hard enough 
with a routinely achievable quality factor of 10*.^^ The 
double dots are defined using electrostatic gates on the 
surface of the heterostructure where the 2DEG is located 
typically 100 nm below the surface^^. 

The resonator is coupled to the DDs magnetically by 
placing a Co micro-magnet on the resonator as depicted 
in fig. [I] with, fioM = 1.8 T, where fiQ and M are the 
permeability of free space and magnetization of Co re- 
spectively. The Co micro-magnet has a length 500 nm. 



width 200 nm and thickness 100 nm and is engineered 
so that the magnetization vector remains parallel to the 
easy axis for small values (~ 50 mT) of the external 
tuning field, Bq. The left and right dots are located 
a distance, c?i = 150 nm and d2 = 250 nm, from the 
mechanical resonator respectively. Under these condi- 
tions, with the nano-mechanical beam at resonance, the 
total field including the external tuning field is of the 
form, B = (Bq + Bzn + SBznSi'n{ujrt))z — BxnSiTi{uJrt)x. 
The stray field of the Co micro-magnet has been calcu- 
lated analogously to the method outlined in ref.^'^. The 
in-plane (Zeeman) magnetic field, Bz, arising from the 
above described magnetic element is shown in fig. [2^. 
At the QD locations, the in-plane magnetic field at the 
left (right) dot is, B^i(2) = 94 mT (36 mT). The in- 
plane component of the magnetic field exhibits almost 
no temporal dependence {5Bzn ^ Bzn) even though the 
field originates from the magnetic element placed on the 
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resonator because the beam oscillation is much smaller 
than the thickness of the micro-magnet. The proximity 
of the magnetic element to the DDs produces a strong 
in-plane field gradient of order, {B^i — 3^2) /{di — ^2) ^ 
0.6 T//im which allows the nano-mechanical oscillator to 
addressably couple to the electron spin in cither quan- 
tum dot as discussed below. The spatial dependence of 
the transverse magnetic field, B™^^, whilst the beam is 
at maximum deflection is shown in fig.[2]3. In the inset of 
fig. [Dd, the time dependence of the transverse magnetic 
field, B^, at the locations of the DDs whilst the beam 
is at resonance is shown. The amplitude of this field at 
the location of the left (right) dot is - 3.0 mT (0.8 mT) 
which is more than sufficient to flip the electron spin in 
either dot via ESR. As can be seen from fig. [5^ and b, 
the parameter, do, is of critical importance as this gov- 
erns the degree to which the DDs couple to the resonator 
via the magnetic element. 



A. Hamiltonian 

For two-electron QDs, the relevant spin states are the 
spin singlet and triplet states of the (1,1) charge state 
(|5'),|To),|T+),|T_)) and the spin singlet of the (0,2) 
charge state (| (0,2)5')) where the label (m,n) refers to 
the number of electrons confined on the left and right 
dots2i. In the presence of finite interdot tunnelling, T, 
the (1,1) and (0,2) singlets are coupled allowing electrons 
to be moved between dots when the detuning parameter, 
e, is varied (by pulsing the QD gate voltages). The en- 
ergy, = 5/iB(S^i - -B^2)/2, couples the \S) and \Tq) 
states whilst |T+) and |T_) remain eigenstates with Zee- 
man energy, ±Ez = ±gnBiBzi + + 2Bo)/2, where, 
g, is the electron g-factor {\g\ = 0.44 for GaAs) and, 
/is, is the Bohr magneton. For the stationary resonator, 
the Hamiltonian describing the two electron states in the 
presence of the above described in-planc in-homogeneous 
magnetic field in the {|(0, 2)5), |5), \To), \T+), |T_)} basis 
is 



-e V2T 

V2T -A^ 

H= -A^ 0. (1) 

-Ez 

Ez 

In fig. [3^, the eigenenergies of H are plotted as the de- 
tuning is swept between the (1,1) and (0,2) charge states 
across the degeneracy point (e = 0). For large positive 
detuning, the ground state is |(0, 2)5). Close to e = 0, T 
and A^ mix the |(0,2)5), |5) and \To) states. For large 
negative e, | (0,2)5) uncouples from |5) and \To) to be- 
come the highest energy state. In the absence of a field 
gradient (A^ = 0), |5) and \To) would also be eigenstates 
with splitting equal to the exchange energy, J « T^/e. 
However, because A^ ^ J, the actual eigenstates are the 
superpositions, ((|5) ± \Tq))/\/2), which correspond to 



the states |ti) (spin up and spin down at the left and 
right QDs respectively) and ||t) respectively. Moreover, 
the field gradient in the regime of large negative detuning 
results in asymmetric separation of jT-f) and |T-_) states 
about I'll) where the energetic separation between the 
Iti) and |T+> = ITT) states is, E^2 - 5Ms(So +5^), the 
Zeeman energy of the right spin. 
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FIG. 3: (a) (colour on-line) The two-electron spin states as 
a function of detuning energy with interdot tunnelling, \/^T 
= 1 /leV with the external tuning field set to zero. At large 
negative detuning, the transverse ESR field {Bx2) resonantly 
couples the and |r+) states (blue line) when the con- 
dition, g^B{Bz2 + Bo) = hur, is satisfied, (b) The pulse se- 
quence consisting of preparation in the |(0, 2)S} (A), resonant 
coupling of the and |r+) states for time, tesr, (B-^C) 
and readout (D) steps. For, tesr ~ Ttt, the final state at 
point (D) is the triplet |r+). 



III. SPIN ISOLATION AND ESR 

The nano-mechanical resonator can be used to ma- 
nipulate the electron spin states in the DDs as follows. 
The resonator is excited at its natural frequency corre- 
sponding to an energy, huJr, the equivalent magnetic field, 
TiLjj./ g^B is designed to be much greater than GaAs hy- 
perfine magnetic field, Bf^ 1 mT. The in-plane tuning 
magnetic field is then used to adjust the Zeeman spin 
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FIG. 4: The expected QPC measurement for the occupation 
of the |1, l)r+ state as a function of the field, Ez2/ffA's, for 
long ESR mixing time, tesr, when the right QD is resonantly 
coupled to the mechanical element. The resolution of the 
ESR peak is reduced and broadened by the increasingly strong 
nuclear field. The range of fields presented here are (-Bjv) = 
0, 1, 5 mT, with the expected hyperfine field in GaAs being 
ImT. The solid lines represent Q = 10* while the dashed lines 
represent Q — 10^, illustrating the inverse dependance of the 
QPC line shape on the quahty factor of the resonator. Here 
the incoherent 7 process is negligibly small. 



splitting to couple the right QD and the mechanical ele- 
ment by setting — 'hur- That is, a tuning field corre- 
sponding to the energetic difference between the Zeeman 
spin splitting at the right QD and the mechanical oscil- 
lator results in resonant coupling. For a 200 MHz res- 
onator with the above values for the in-plane magnetic 
field, the tuning field will be -3.5 mT. With the mechan- 
ical element at resonance, the 1(0,2)5') state is prepared 
(A) which is then separated into the |ti) state by the in- 
plane magnetic field gradient between the two QDs and 
by using rapid adiabatic passage to change e from a pos- 
itive to a negative value (A^B)^^. With the right QD 
resonantly coupled to the resonator, the transverse ESR 
magnetic field couples only the |ti) and |T+) states at 
a rate proportional to Yix2 = gi^BBx2- The left spin is 
unaffected by the ESR field due to the large difference 
in Zeeman field between the QDs. Assuming an infinite 
quality factor, the Hamiltonian in these conditions is 



H 



ESR 



-^zi<7z2l^ + Sj;2(wr, ^d, (3)sin(a;dt)cra;2/2, (2) 



where cr„2 — (px2^ (^y2,o'z2) are the Pauli matrices of the 
spin located in the right QD, uj^ is the natural frequency 
of the resonator, and ujd is the driving frequency of the 
resonator. The magnitude of the Tix2 field changes as a 
function of the quality factor, and w^, according to its 
proportionality to the displacement of the resonator 
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v/(a;2_^.)2+4^2^./(Q2)' 



(3) 
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where C is a proportionality constant. For later use we 
define 6uj = uj?. — ui'i . 



FIG. 5: A contour plot of the expected QPC measurement 
for the occupation of the |1, l)r+ state as a function of the 
field, Ez2/<? /is, and the driving frequency of the resonator cOd, 
for long ESR mixing time, tesr, and hyperfine interaction 
(Bjv) = ImT. The inset shows the same response for the 
zero hyperfine interaction. The horizontal arm of the 'star' is 
defined by Tiz2/'h — u>d ~ 0, as illustrated by the superimposed 
line in the inset. When the hyperfine interaction is on, the 
horizontal arm is broadened considerably, with only a small 
exaggerated response at the ujr = uJd point. Here, white is 
maximum (0.42 for the main figure, ~ 0.5 for the inset), black 
is minimum. 



The coupled system is held in this state for time, tesr, 
until steady state has been achieved (B— >C). After ma- 
nipulation, rapid adiabatic passage is used to change e 
from negative to positive values for charge sensing mea- 
surement via the quantum point contact, QPC (C^D). 
The |(1,1)T_|_) state remains in the spin-blocked config- 
uration whereas the |ti) tunnels directly to |(0, 2)5). A 
schematic of this pulse sequence is shown in fig. Af- 
ter manipulation, the spin state is a mixture of |(0, 2)S) 
and 1(1, 1)T-|_) states (at point D). In this condition, the 
QPC current, Iqpc, will be a mixture of the correspond- 
ing (0,2) and (1,1) charge states. At off resonance, the 
spin configuration at (D) will be |(0, 2)5) with (0,2) QPC 
signal. Therefore, by sweeping the tuning field across res- 
onance, a peak in the Iqpc is expected. 



IV. MASTER EQUATION OF ESR 
HAMILTONIAN 



The response of the QPC can be modeled using a mas- 
ter equation model for Hesr- To robustly model the 
coupled system, the various mechanisms by which the 
electron spins can decohere are considered. The |(0, 2)5) 
is initialized at point A in the schematic of the pulse se- 
quence shown in fig. ^p. On experimental time scales. 
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this state is highly robust and no spin dephasing takes 
place. Concurrently the mechanical element is at reso- 
nance and thermal fluctuations in its motion can result 
in the ESR field smearing out i.e. the time dependence 
of the ESR field shown in the inset of fig. [2}d will result 
in small amplitude oscillations being super-imposed on 
the sinusoidal ESR response. However, the amplitude of 
these thermal fluctuations can readily be extracted from 
the equipartition theorem and for the operating tempera- 
ture of this experiment (~50 mK) the thermal amplitude 
is expected to be '^O.l pm which is more than 4 orders of 
magnitude less than the amplitude of the beam. There- 
fore, spin dephasing due to thermal noise in the mechani- 
cal element is negligible. Furthermore, the Curie temper- 
ature of the Co micro- magnet is 1394 K which is much 
greater than the operating temperature of the proposed 
experiment. Hence spin dephasing due to fluctuations in 
the ferromagnetic properties of Co are again negligible. 
In general electrical noise (the detuning parameter, e, is 
varied by pulsing the voltages of the QD gates) induced 
fluctuations in the exchange energy, J (via A and T 
noise) lead to extra mixing between the |(0, 2)5), \S) and 
\Tq) states. However, because the spins are manipulated 
at large detuning, gate voltage noise can be neglected be- 
cause dJ/ dA ~0 in this regime, i.e. the states are to first 
order insensitive to electrical noise. As demonstrated in 
refi^, spin-orbit coupling, which is an extra decoherence 
channel can be neglected in GaAs quantum dot system 



for the proposed experimental timescales. The remain- 
ing main decoherence mechanism for the electron spin is 
the Hyperfine interaction with nuclear spins within the 
quantum dots^^. We ignore the hyperfine interaction for 
the electron in the left dot by assuming the mean nuclear 
magnetic field is much smaller than the Zeeman splitting, 
thus the spin state is not flipped. In addition, because of 
the large energy difference between the two dots, we can 
ignore inelastic singlet-singlet transitions. 

Under these conditions, the master equation describing 
the spin states of the electron in the right dot is. 



= -i[HESR + HN.p[t)]+L[p{t% (4) 

where Hf^ — —g.fiB{BN-<^2) is the hyperfine interac- 
tion and L[/5(i)] describes the incoherent Markovian spin 
relaxation at rate 7. We reach the steady state solution 
by taking t — > cxd, ^^^^ = 0, and thus inverting to find 
p{t ^ 00) = [-ifi[HESR + HM,p{t)] + Assuming 
the rotating wave approximation, we are able to obtain 
an analytical solution for the occupation of the |(1, 1)T+) 
and 1(0, 2)5") states. 

For static hyperfine nuclear fields Bm, the steady state 
of the density matrix elements for the diagonal |(1, l)T-f), 
1(0,2)5'), and the off-diagonal ((1, l)T+|p|(0, 2)5), states 
are (omitting the notation t — )■ cxj). 



P(1,1)T+ 
P(0,2)S 
(0,2)5 



{4.Bx,N + S,2) 



32(Bi^ + 

1 - P(1,1)T+ 



52 



X,N 



iBY,N) + S^2)(4Sz,Ar + 2hAuj - ihj) 



l6{Bj^ + By jv + 2i?2 n) + 16Sj;2-Bx^ 



N 



r 



(5) 
(6) 
(7) 



where Au = T,z2/^ — ojd is the Rabi frequency, 7 is a 
phenomenological damping of the spin state, and Bi^jq 
is the nuclear field in the ith direction. The off-diagonal 
term is in the rotating frame. We assume that the nu- 
clear field is static on the h/T,x2 and I/71 time scales, 
but is fast on the time scale of collating a clear QPC 
measurement. Thus, in the data presented in fig. 31 we 
integrate over a normal distribution of possible nuclear 
spin orientations^ with mean value (-Bat) = \/3(-Bi,Ar), 
and a standard deviation of \/2{B]y), so 



1)T+/Bis 



p(i,i)T+ n 



(2(B,,jv)0F) 



-dB, 



N- 



(8) 



Consequently, (B^) is the important quantity in deter- 



mining the affect of the nuclear field on the ESR measure- 
ment. In fig. m the effects of the hyperfine coupling on 
the QPC measurement are shown. For (Bn) < 'Ex2 the 
on resonant response of the QPC is clearly visible. For 
{Bq^n ^ Sa;2) the on resonance peak becomes smeared 
out and eventually obliterated. 



V. ESR LINE WIDTH AND RESONATOR 
QUALITY FACTOR 

Ignoring the hyperfine interaction and the stochastic 
decoherence rate of the spin (7 = 0), we can see very 
clearly how the line width of the ESR response depends 
on the Quality factor of the resonator. The ESR response 
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FIG. 6: Contour plots of the expected QPC measurement for 
the occupation of the |1, 1}T+ state as a function of the field, 
Tjz2/gfJ'B, and the driving frequency of the resonator ua, for 
zero hyperfine interaction, and a range of Q factors. Changing 
the Q factor alters the range of the vertical arm {ujd — uJr~0) 
of the star. In the top two figures, white is maximum. In the 
bottom two figures, black is maximum (for clarity). 



P(1,1)T+ 



2(l + §^Ac<j2)' 



(9) 



Hovifever, as mentioned earlier, the field ^,^2 is both Q 
and uJd dependant. This gives the following Lorentzian 
form with two maxima. 



P(1,1)T+ 



1 



2(1 



C2 



(10) 



Thus the line width of the ESR response, obtained by 
sweeping the tuning field, is 
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ESR 



QC 



(11) 



where C again is the constant defined by properties of 
the resonator (see Eq. [3] ) . When the resonator is driven 
at resonance = w^, the ESR line width is 7^ ~ Q, the 
inverse of the resonators own line width, Eq. [3] . This 
can be understood by the fact that the reduced motion of 
the resonator, because of a lower Quality factor, implies 
a reduced maximum T,x2 field, which in turn reduces the 
steady state probability for the spin to be in the | (1, 1)T+) 
state when the Rabi tuning Aoj is off-resonance. Fig. [6] 
illustrates the effect of changing Q on the QPC measure- 
ment. 



In addition, for large detuning of the resonator driv- 
ing frequency from the natural frequency Suj = uj'^ — cu'^ , 
we find that 7/^ w 1/ y/^uj'^ ~ ^r)- Hence the Alj peak 
becomes very sharply defined and almost independent 
of the quality factor. Thus the transduction of the res- 
onator motion by ESR is advantageous because it allows 
us to read out both the quality factor Q and intrinsic 
properties of the resonator via the constant C . Further- 
more, for large detuning uj"^ ^ ^"i, the spin rotation can 
be performed with a very narrow line width. 

Alternatively, by fixing the tuning field at Bq — 
hujr/giiB — Bz2 {Bz2 can be measured in principle by 
DD spectroscopy^^) and sweeping the frequency of the 
mechanical oscillator about the resonance point, 5lo — 
Aw = 0, a peak in Iqpc is again seen. But now the cor- 
responding line shape is the convoluted response of the 
nano-mechanical oscillator and the two-electron system. 
If the tuning field is is moved away from the LUr/gfJ^B 
point, both maxima in the double Lorentzian are visible. 
In this case, the two peaks both have a cod dependant line 
shape, which dominates over the Q dependance. 

Additionally, moving beyond the steady state, the me- 
chanical resonator provides a local ESR field which en- 
ables the electron spin states to be manipulated. For 
example, by modulating the waiting time (B^C) to pro- 
vide a TT pulse so that, tesr = = '2TTh/'Ex2 — 0.2 
fis , then the final state at readout is the triplet 
By employing a larger tuning field, the states |ti) and 
|r_) can also be resonantly coupled so that the spin of 
the left dot is fiipped given the tuning field is less than 
the Co switching field so that the magnetisation does not 
reverse. 

In the above described coupled system, actuation of 
the mechanical element is achieved via the piezoelectric 
effect between the top gate defined over one of the clamp- 
ing points of the beam and the shallow buried 2DEG2-. 
The mechanical resonance is then excited by applying an 
AC voltage to the top gate at the mechanical resonance 
frequency. Crucial to the success of this scheme is the 
in-plane alignment of the micro-magnet and the shallow 
buried 2DEG on the mesa on which the DQDs are de- 
fined. This can be achieved by applying a DC voltage 
to the top gate on the mechanical element which will in- 
duce a static deflection of the mechanical beam enabling 
alignment to be achieved as visualised in the inset of 
fig. [2^. Simultaneously a small AC modulation will drive 
the beam to mechanical resonance. 



VI. CONCLUSIONS 

A two-level system realized in a two-electron double 
quantum dot is a highly sensitive quantum detector. 
Magnetically coupling the double quantum dots to a high 
frequency nano-mechanical resonator will allow the dis- 
placements of the resonator to be detected. Tuning the 
Larmor frequency to match the resonators natural fre- 
quency, the electron spin in one of the dots can be flipped. 
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By simultaneously measuring the charge state of the dou- 
ble quantum dots, this transition can be detected thus en- 
abling the natural frequency of the mechanical resonator 
to be determined. 

The authors NL, IM and MP-L have contributed 
equally to this manuscript. The authors also acknowl- 
edge A. Blais for useful discussions and for critically 
reading the manuscript. S. Camou, S. Etaki, C. Buizert 



and H. Inagaki are thanked for their help and encour- 
agement. We acknowledge financial support from the 
Grant-in- Aid for Scientific Research A (No. 40302799), 
the Special Coordination Funds for Promoting Science 
and Technology, MEXT, SORST-JST, DARPA grant no. 
DAAD19-01-1-0659 of the QuIST program and JSPS 
KAKENHI(16206003) and JSPS grant 17-05761. 



* Electronic address: llamb ert @asone . c . u-t okyo . ac. j p | 
^ H. G. Craighead, Science 290, 1532 (2000). 

^ A. Cho, Science 299, 36 (2003). 

3 M. Roukes, Phys. World 14, 25 (2001). 

■* R. G. Knobel and A. N. Cleland, Nature 424, 291 (2003). 

^ M. D. LaHaye, O. Buu, B. Camarota, and K. C. Schwab, 

Science 304, 74 (2004). 
^ A. Naik, O. Buu, M. D. LaHaye, A. D. Armour, A. A. 

Clerk, M. P. Blencowe, and K. C. Schwab, Nature 443, 

193 (2006). 

M. R. Geller and A. N. Cleland, Phys. Rev. A 71, 032311 
(2005). 

* F. Xue, Y. D. Wang, C. P. Sun, H. Yamaguchi, and 
K. Semba, arXiv:cond-mat/0607180vl. 

^ S. M. Carr, W. E. Lawrence, and M. N. Wybourne, Phys. 

Rev. B 64, 220101 (2001). 
^° S. Savel'ev, X. Hu, and F. Nori, arXiv:cond-mat/0412521. 
" A. Wallraff, D. Schuster, A. Blais, L. Frunzioand, R. S. 

Huang, J. Majer, S. Kumar, S. M. Girvin, and R. J. 

Schoelkopf, Nature 431, 162 (2004). 
12 R. C. Ashoori, Nature 379, 413 (1996). 
1^ S. Tarucha, D. G. Austing, T. Honda, R. J. van der Hage, 

and L. P. Kouwenhoven, Phys. Rev. Lett. 77, 3613 (1996). 
1* M. Ciorga, A. S. Sachrajda, P. Hawrylak, C. Gould, P. Za- 

wadzki, S. Jullian, Y. Feng, and Z. Wasilewski, Phys. Rev. 

B 61, R16315 (2000). 
1^ D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 

(1998). 

J. Petta, A. Johnson, J. Taylor, E. Laird, A. Yacoby, 
M. Lukin, C. Marcus, M. P. Hanson, and A. Gossard, Sci- 
ence 309, 2180 (2005). 
1^ F. H. L. Koppens, C. Buizert, K. J. Tielrooij, I. T. Vink, 



K. C. Nowack, T. Meunier, L. P. Kouwenhoven, and 
L. M. K. Vandersypen, Nature 442, 766 (2006). 
1* F. Xue, L. Zhong, Y. Li, and C. P. Sun, arXiv:quant- 
ph/0602219v2. 

1^ D. Rugar, R. Budakian, H. J. Mamin, and B. W. Chui, 
Nature 430, 329 (2004). 

2° J. M. Elzerman, R. Hanson, L. H. W. van Beveren, 
B. Witkamp, L. M. K. Vandersypen, and L. P. Kouwen- 
hoven, Nature 430, 431 (2004). 

21 At large oscillation amplitudes bi-stability occurs where 
the effective spring constant of the beam increases and the 
resonance peak migrates to higher frequencies simultane- 
ously developing asymmetry at the expense of quality fac- 
tor. 

22 J. M. Elzerman, R. Hanson, J. S. Greidanus, L. H. W. 
van Beveren, S. D. Franceschi, L. M. K. Vandersypen, 
S. Tarucha, and L. P. Kouwenhoven, Phy. Rev. B 67, 
161308 (2003). 

2^ M. Pioro-Ladriere, Y. Tokura, T. Obata, T. Kubo, and 
S. Tarucha, Appl. Phys. Lett. 90, 024105 (2007). 

24 J. M. Taylor, H.-A. Engel, W. Ur, A. Yacoby, C. M. Mar- 
cus, P. ZoUer, and M. Lukin, Nature Physics 1, 177 (2005). 

25 A. V. Khaetskii and Y. V. Nazarov, Phys. Rev. B 64, 
125316 (2001). 

2® W. A. Coish and D. Loss, Phys. Rev. B 72, 125337 (2005). 

2^ F. H. L. Koppens, C. Buizert, L T. Vink, K. C. Nowack, 
T. Meunier, L. P. Kouwenhoven, and L. M. K. Vander- 
sypen, J. Appl. Phys. 101, 081706 (2007). 

2® S. D. Barrett and T. M. Stace, Phys. Rev. Lett. 96, 017405 
(2006). 

2^ I. Mahboob, S. Miyashita, and H. Yamaguchi, unpub- 
lished. 



